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ABSTRACT 

We have studied the filaments extracted from the column density maps of the 
nearby Lupus 1, 3, and 4 molecular clouds, derived from photometric maps observed 
with the Herschel satellite. Filaments in the Lupus clouds have quite low column 
densities, with a median value of ~1.5xl0^^ cm“^ and most have masses per unit 
length lower than the maximum critical value for radial gravitational collapse. Indeed, 
no evidence of filament contraction has been seen in the gas kinematics. We find 
that some filaments, that on average are thermally subcritical, contain dense cores 
that may eventually form stars. This is an indication that in the low column density 
regime, the critical condition for the formation of stars may be reached only locally 
and this condition is not a global property of the filament. Finally, in Lupus we find 
multiple observational evidences of the key role that the magnetic field plays in forming 
filaments, and determining their confinement and dynamical evolution. 

Key words: ISM: clouds - stars: formation. 


1 INTRODUCTION 

Filaments are ubiquitous structures present in both low- 


and in high-mass star-forming regions (e.g., Hatchell et al. 
2005| [Goldsmith et al.||2008| [Andre et al.||2010| [Molinari 


et al.pOlO l. They were recognised more than 30 years ago 
I Schneider & Elmegreen 19791 but only recently the Her- 


sc/ie0Space Observatory has shown the omnipresence of the 
filamentary structures throughout star-forming regions and 
their strict association with the dense cores from which new 
stars will form (e.g., Andre et al.|2010 Molinari et al.|2010 


Polychroni et al. 20131. The Herschel results have focused 
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special attention on the formation mechanism of filaments, 
and in understanding their role in the star formation process, 
with several studies either analysing interesting filaments, 
or adopting a more statistical approach. These studies, al¬ 
though far from answering the key questions, have started 
to reveal important properties. For example, filaments in 
molecular clouds are characterised by having a narrow dis¬ 
tribution of widths, with a median value of 0.10±0.03 pc 
( Arzoumanian et al.pOlT l. This corresponds, within a fac¬ 
tor of 2, to the sonic scale below which interstellar turbu¬ 
lence becomes subsonic in diffuse gas. This similarity in scale 
supports the argument that filaments may form as a result 
of the dissipation of large-scale turbulence. Lower column 
density filaments are thermally subcritical and gravitation¬ 
ally unbound, and are probably confined by external pres¬ 
sure, while higher column density filaments are thermally 
supercritical and gravitationally bound and show evidences 
of gravitational collapse (e.g., Schlsano et al. 2014[ Arzouma- 


nian et al.|2013 Nguyen-Luong et al.|201^Hill et al.|2012 |. 
A kinematical study conducted on the B211/213 filament 
in Taurus ( Hacar et al.||2013 l, revealed that some filaments 
are actually composed of several distinct, coherent velocity 
components of gas. In the more massive filaments, spectral 
observations indicate a more dynamical condition of the gas, 
revealing velocity gradients indicative of collapse and con¬ 
traction (Schneider et al. 112010 Kirk et al. 11201 3 Peretto et 

[iT][^ -- -- - 

The Lupus complex is an interesting region in which 
to study filaments. It is one of the nearest low-mass star¬ 
forming regions (distance between 150 pc and 200 pc. 


Comeron 2008) and a few long and well defined filaments 


have been identified in it (Rygl et al. 2013 Benedettini et 


al. 2012 Poidevin et al. 2014 Matthews et al. 20141. The 


complex is located in the Scorpius-Centaurus OB associa¬ 
tion, whose massive stars are likely to have played a signifi¬ 
cant role in shaping the diffuse material and in the general 
evolution of the complex. 

In this paper, we present an analysis of filaments ex¬ 
tracted from the column density maps of the Lupus 1, 3, 
and 4 star-forming regions, produced from the photometric 
observations carried out with the Herschel satellite. In Sect. 
2 we describe the observations and the data reduction. The 
analysis of the filament properties and of the probability dis¬ 
tribution function is presented in Sect. 3 ad 4, respectively. 
The results are discussed in Sect. 5 and the main conclusions 
are reported in Sect. 6. 


2 OBSERVATIONS AND DATA REDUCTION 
2.1 Herschel photometric maps 

As part of the Herschel Gould Belt Survey (HGB^ [Andre 
et al.|20l0 1 the three subregions of the Lupus complex, Lu¬ 
pus 1,3, and 4, have been observed in five photometric bands 
between 70 pm and 500 pm with the Photodetector Array 
Camera and Spectrometer (PACS; [Poglitsch et al.|2010 ) and 
the Spectral and Photometric Imaging Receiver (SPIRE; 
Griffin et al. 2010| o n-board the Herschel Space Observa¬ 
tory (Pilbratt et al. 2010). The observations were carried 


out in the parallel observing mode with a scanning velocity 
of 60 arcsec/sec and each zone was observed twice with or¬ 
thogonal scanning directions. The spatial resolution in this 
observing mode is 7.6” at 70 pm, 11.5” at 160 pm, 18^2 at 
250pm, 24”9 at 350 pm and 36”3 at 500 pm. A first ver¬ 


sion of the Lupus maps was published by Rygl et al. (2013) 
who used the romagal pipeline ( Traficante et al.|2011 |. In 
this paper we present a new reduction of the maps obtained 
with a different pipeline and updated calibration files. The 
pipeline for the production of the PACS maps (at 70 pm and 
160 pm) consisted of the following steps: the Level 1 data 
products were generated with Standard Product Generation 
(SPG) version 10 within the Herschel Interactive Processing 
Environment (hipe[^; Level 1 data were arranged in time or¬ 
dered data (TOD) and exported in fits format; the TODs 
of the two orthogonal observations were ingested into the 
UNIMAP map maker ( Piazzo et al.|[2015 l to produce the fi¬ 
nal maps. The SPIRE data (at 250 pm, 350 pm, and 500 
pm) were reduced using modified pipeline scripts of hipe 
version 10. Data taken during the turnaround of the satel¬ 
lite were included in order to assure a larger area covered by 
the observations and enough overlap. Because the observed 
Galactic regions are dominated by extended emission of the 
interstellar medium, extended emission relative gain factors 
were applied to the bolometer time-lines. The resulting Level 
1 contexts for each scan direction were then combined using 
the naive map maker in the destriper module. 

The final absolute flux calibration for both PACS and 
SPIRE maps was performed by adding a zero-level offset 
derived by comparing the Herschel data with the Planck 
and IRAS data of the same area of the sky. A dust model 
has been assumed for extrapolating the flux at the Herschel 
wavelengths ( Bernard et al.j 20101. The final flux accuracy 
is better than 20 per cent. Notably, the Lupus 1 maps were 
affected by stray Moonlight, visible as a bright vertical band 
in each image. We removed this contamination by evaluating 
the stray-light contribution from the difference between the 
observed map (affected by stray-light) and the model used 
for the flux calibration. This procedure worked well but not 
for the 70 pm image for which the correction could not be 
applied. The 70 pm map, however, was not used for the 


production of the column density maps (see Sect. 2.31, that 
therefore are corrected from the stray-light effect. 


2.2 CS (2—1) line maps 

A CS (2-1) line map of the Lupus 1 and 3 clouds was 
acquired at the 22-m Mopra radio telescope using the On 
The Fly observing mode with the narrow band mode of the 
UNSW-Mopra Spectrometer (UNSW-MOPS) digital filter- 
bank back-end. The velocity resolution at 97.981 GHz is 
0.1 km s“^ the typical rms of the single spectrum is ~ 0.4 
K and the Half Power Beam Width is 35 arcsec. Part of 
these maps have been already published in [Benedettini et| 
al. (2012), where details regarding the data reduction can 


http: / / gouldbelt-herschel.cea.fr 


be found. In May 2012, we acquired new CS (2-1) spectra 


^ HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, The NASA Herschel 
Science Center, and the HIFI, PACS and SPIRE consortia 
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Table 1. Parameters of the Lupus clouds. 



Lupus 1 

Lupus 3 

Lupus 4 

Distance (pc) 

150 

200 

150 

area above Av >2 mag (deg^) 

0.32 

0.19 

0.22 

Mass above Ay >2 mag (Mq) 

163 

162 

106 


toward the Lupus 3 cloud extending the coverage of the re¬ 
gion with respect to the previous published map. The CS 
(2-1) maps cover only the brighter and denser part of the 
Lupus 1 and 3 clouds, approximately the region above the 
4 mag of visual extinction. 


2.3 Column density maps 


The column density maps of the three Lupus clouds are 
shown in Fig.[^[^and[^ The maps were produced by apply¬ 
ing a grey-body fit in each pixel to the Spectral Energy Dis¬ 
tribution (SED) of the four surface brightness values at 160 
/rm, 250 /rm, 350 /rm and 500 ^m smoothed to a common 
spatial resolution of 36 arcsec, i.e. the resolution of the 500 
pm map and resampled with a pixel size of 14 arcsec. The 
spatial resolution of the column density maps corresponds to 
a linear scale of 0.026 pc for Lupus 1 and 4 (distance of 150 
pc) and 0.035 pc for Lupus 3 (distance of 200 pc). For the 
grey-body htting, we assumed a dust opacity of K25o=0.144 
cm^g“^ (already accounting for a gas-to-dust ratio of 100), 
a grain emissivity parameter /3=2 (cf. Hildebrand 1983) and 
a mean molecular weight p=2.8, values adopted as the stan¬ 
dard by the HGBS consortium. Previous works (e.g., 1^ 
et al.|20T4 Sadavoy et al.|2013 1 estimated that the column 
densities derived with these assumptions are accurate to 50 
per cent for dense cores. The H 2 column density towards the 
Lupus clouds is typically low, with the denser region corre¬ 
sponding to visual extinction between 3 mag and 6 mag and 
a very small portion of the clouds at > 7 mag. 

We estimated the total mass of the clouds by summing 
the mass inside the Av=2 mag contour, using the conversion 


formula AiH 2 = 9.4x10^“ Av (Bohlin et al. 19781 and the 


derived values are listed in Table [Ij We obtained masses 
lower than the estimates previously quoted by |Rygl et al.| 
(20131 from the Herschel maps, namely 830 Mq (the map 
was not corrected for the stray-light), 570 Mq and 500 Mq 
for Lupus 1, 3 and 4, respectively. The discrepancy is mainly 
due to the fact that they integrated over an area larger than 
our area. In fact, the original five bands maps from which 


Rygl et al. (20131 derived their column density maps were 
produced with a different map-maker and older calibration 
files and they used a different opacity law. 

In the literature, we have found a variety of estimates 
for the mass of the Lupus clouds, derived by using extinction 
or molecular line maps, spanning a large range of values. A 
direct comparison between our mass estimates and the lit¬ 
erature values is not so trivial since they are derived from 
maps covering different area of the clouds; moreover the star 
counts extinction maps and the CO maps have spatial reso¬ 
lution of a few arcminutes and are usually spatially under¬ 
sampled. Generally speaking, our masses are more similar to 


the values derived from ^CO (2-1) line maps (Tothill et al. 


20091 and are lower than previous estimates from other ex¬ 
tinction maps (Merin et al.|200^ Cambresy|1999 1. In order 


to do a more accurate comparison we produced new visual 
extinction maps by using the star count method of |Cambresy| 
(19991 applied to near-infrared (NIR) data. These new ex¬ 
tinction maps cover the same region mapped with Herschel 
and have a spatial resolution of 2 arcmin. In Fig. |^we com¬ 
pare the visual extinction derived from the NIR data with 
that derived from Herschel data, in all the common spatial 
pixels, smoothed at the common resolution of 2 arcmin. We 
find a good correlation between the two values, but above 
the threshold of ~ 2 mag the Herschel values are lower of 
a factor of ~ 2 than the star counts values, although the 
higher spatial resolution of the original Herschel maps. Sim¬ 
ilarly, column density higher than the Herschel values have 
been derived from Planck data at the even coarse spatial 
resolution of 10 arcmin (Planck Collaboration 20151. It is 
worth noting that the Herschel maps detect the dust emis¬ 
sion and the extinction maps derived from these measure¬ 
ments require the assumption of the dust opacity, the grain 
emissivity and the dust to gas ratio. It is therefore possi¬ 
ble that the observed discrepancy arises from a wrong as¬ 
sumption in the dust parameters that we used in deriving 
the column density maps. Values of visual extinction more 
similar to the star counts values could be reproduced from 
the Herschel data halving the opacity at the reference fre¬ 
quency, i.e. assuming K25 o= 0.072 cm^g~^. In fact, several 
studies demonstrated that the dust opacity is not constant 
in the diffuse ISM but it is correlated to the Hydrogen col¬ 
umn density and more in general to the volume density. For 
example, Planck Collaboration (2011aI found K 250 = (0.026- 
0.069) cm^g“^ for high latitude diffuse ISM where the gas is 
mainly in atomic form, while in the Taurus molecular cloud 
Planck Collaboration (2011b I found K25 o= 0.049 cm^g“^ for 


the atomic phase gas and K25 o= 0.099 cm^g~^ for the molec¬ 
ular gas. Martin et al. (2012 1 found K25 o=( 0.086 - 0.17) 
cm^g”'^ for a sample of 14 regions near the Galactic plane 
toward the Vela molecular cloud, mostly selected to avoid 
regions of high column density (N(H) > 10^^ cm“^). Fi¬ 
nally, Roy et al. (20131 analysing the Herschel data of the 
Orion A molecular cloud found that the dust opacity in the 
region increases from 0.043 to 0.129 cm^g“^ as function of 
jY(H)0'28 Qygj. range of A'(H) from 1.5 xlO^^ cm“^ to 
50 X 10^^ cm~^. Looking at these results, the dust opac¬ 
ity of K250=0.144 cm^g“^ assumed in the HGBS consortium 
for the SED htting of the dense cores could be too high for 
the calculation of the column density map in the Lupus com¬ 
plex, which column density is typically lower than the values 
observed towards dense cores. However, for consistency with 
the other studies of the HGBS project, in the following anal¬ 
ysis we assumed the HGBS standard, bearing in mind that 
the derived column density values could be underestimated 
by a factor of two. 


2.4 Filament detection 

As Herschel has largely shown, molecular clouds are mainly 
structured in hlaments and we see such structures also in 
the Lupus clouds. In order to characterise the hlamentary 
structure of the Lupus clouds, we applied an algorithm that 
uses the eigenvalues of the Hessian matrix of the column 
density map to identify the regions belonging to hlaments. 
At the hrst order, the border of the Hlament is identihed 
as the place with the maximum variation of the contrast of 
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15 40 "' 00 ^ 35 ^^ 00 ^ 

R.A. (J2000) 

Figure 1. Column density map at 36 arcmin resolution derived form Herschel data of the Lupus 1 cloud. The black contours are Av = 
1 mag and 2 mag, respectively. The white lines identify the spines of all the branches of the filaments. 


the column density field. Doing this, the pixels selected will 
always belong to regions where the curvature of the bright¬ 
ness profile in the maximum curvature direction is within 
the convexity region. This implies in general a conserva¬ 
tive identification of the filament region, because it would 
approximate the contribution of the wings of the filaments 
since we do not include the region where the emission profile 
changes concavity before joining the background emission. 
In order to recover, at least part of, the external wings of 
the filaments we applied a morphological dilation operator at 
the first order border identified with the thresholding of the 
curvature matrix. The spine of the filament is then derived 
by applying a morphological operator to thin the pixels in¬ 
side the borders, connected with a Minimum Spanning Tree 
(MST) technique. Because of the interwined nature of the 
filamentary structures, each filamentary region is composed 
by one or multiple spines connected in nodal points; we call 
branch a single portion of a filament connecting two nodes. 
A detailed description of the filament detection algorithm 


and its performance can be found in Schisano et al. (20141. 


The filamentary structures identified by the algorithm 


depend on the adopted threshold on the eigenvalue map. The 
eigenvalues are a measure of the region contrast with respect 
to its surrounding, hence the selected threshold determines 
the lowest contrast that an elongated region should exhibt 
to be considered belonging to a filamentary structure. Low 
threshold values trace low contrast regions, including also 
several faint structures. We visually inspected the output 
with different thresholds and choose the value equal to 2.8 
time the local standard deviation of the minimum eigenvalue 
map. In fact such a choise allows to identify the most bright 
structures and many faint filament branches departing from 
the bright ones (often called striation). An example is re¬ 
ported in Fig. 1^ where the detected filament spines are ove- 
layed to a portion of the 250 fim map of Lupus 1 cloud. As in 
any thresholding algorithm, we could miss some faint struc¬ 
tures. To quantify such effect we measured the branch con¬ 
trast of the identified filamentary structure defined as the ra¬ 
tio between the average column density along the spine and 
the average column density at the border. The histogram of 
the branches contrast, reported in the nested panel of Fig. 
shows an increasing numbers of branches for decreasing 
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Figure 2. As Fig. for Lupus 3. 


contrasts up to its peak at contrast around 1.2, testifying 
that we recover well also the low contrast part of hlaments. 
Furthermore, we report in Fig. I^the branch contrast as func¬ 
tion of the average spine column density showing that the 
denser filaments tend to have the higher contrast but that 
well defined low density filaments also exist, although our 
census could to be not complete for branches with very low 
column density contrast. 

Once hlaments have been identihed, we calculated their 
physical properties. The width of a hlament, deconvolved 
from the beam, was derived from the median value of the 
full widths at half maximum (FWHM) of Gaussian ht to the 
inner part of the column density prohle, in the direction per¬ 
pendicular to the local spine along all the pixels of the spine. 
The average column density of a hlament was calculated by 
averaging the column density of all the pixels inside the hl¬ 
ament border and the mass per unit length of each hlament 
was obtained by multiplying the average column density for 
the width. We stress that our dehnition entails global prop¬ 
erties for the whole hlamentary structure. However, given 
the complex structure of the longer hlaments, we also cal¬ 
culated the same physical parameters at a smaller scale, i.e. 


for all the branches composing the hlaments. In particular, 
we estimated the direction of the branches from the linear 
interpolation of the points along the branch central spine. 

An example of the radial column density prohle is shown 
in Fig. it is composed by a hat, dense inner region and a 
power-law behaviour at large radii. This form is typical of 
hlaments detected in the Herschel maps of low mass molecu¬ 
lar clouds ( Arzoumanian et al.|2011 Palmeirim et al.|2013 1. 
Analytically, it is described by a Plummer-like function of 
the form 


Af(H2) = 


Nc(H2)Rflat 

[l + (r/Rflat)2]'^ 


( 1 ) 


that derives from an idealised model of a cylindrical hlament 
laid on the plane of the sky with radial density of the form 


p(r) = 


[1 + 


( 2 ) 


where pc is the central density, A^c(H 2 ) is the central column 
density, Rfiat is the radius of the hat inner region, p is the 
power-law exponent at large radii. For the hlament shown 
in Fig. [^the best ht model has Rfiat ~ 0.03 pc and p — 
1.8, fully in agreement with values found in other hlaments 
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Figure 3. As Fig. ^ for Lupus 4. 





Figure 4. Comparison between the visual extinction derived from the Herschel maps and the star counts maps derived from NIR data 
following the method of |Cambrfey| |T999| |. The continuous line identifies the correspondence between the two values, the dashed line is 
the linear fit to the data. The Herschel based values are systematically lower than the values derived from star counts. 


( [Arzoumanian et al.||201l| jPalmeirim et al1]|2013P and sig¬ 
nificantly shallower than the steep p = 4 profile expected 
for unmagnetized isothermal filaments in hydrostatic equi¬ 
librium ( |Ostriker|1964| ). In Fig.we also report the average 
column density along the spine that some authors use as 
indicative of the filament column density and the average 
column density as calculated by us. The two values differ of 
a factor of about 1.5 well below the uncertainty associated 
to the filament column density (see Sect. 3.11. 


3 FILAMENT ANALYSIS 

3.1 Physical properties of the filaments 

In Figs. andwe show the spines of all the branches of 
the detected filaments on top of the column density maps of 
the three Lupus clouds. Lupus 1 appears itself as a long (~ 

5.2 pc) filamentary structure extending from north-west to 
south-east, broken into a sequence of shorter hlaments with 
lengths from 0.1 pc to 2 pc. South-West of the main fila¬ 
mentary structure, a smaller region of dense material exists 
with a roundish form where other filaments have been iden¬ 
tified. Similar to Lupus 1, the Lupus 3 cloud is composed 
of two long hlamentary structures, running almost parallel 
in the east-west direction and fragmented into shorter fila- 
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Figure 5. A zoom into the densest region of the Lupus 1 cloud 
that shows the 250 fira emission (background), the spine of the 
identified filaments (white lines) and the CS (2-1) integrated in¬ 
tensity (cyan contours). Note that the CS (2-1) map does not 
cover the Northern (declination > -33°50') and the Western (R.A. 
< 15^42"*) part of the map. 



10 "“ 10 "' 10 "" 

Average Branch Spine Column Density (cm'^ 


10 "' 


Figure 6. Branch contrast, defined as the ratio of the average 
column density along the spine of the branch and the average 
column density at the border, vs the average spine column density. 
The up-left panel show the histogram of the branch contrast. 


merits. By contrast, the distribution of filaments in Lupus 4 
appears to be more chaotic, without a clearly defined several 
pc filamentary structure, even though several long (up to 3 
pc) filaments running east-west have been identified. 

In Fig. we present the distributions of the physical 
properties of the filaments identified in each of the three 



radius (pc) 


Figure 7. Mean radial column density profile (crosses) of the 
long filament displayed in the Northern part of Fig.[^ The error 
bars show the Icr dispersion of the distribution of radial profiles 
along the brightest part of the filament. The red continuum curve 
shows the best-fit Plummer model described by Eq.^withp = 1.8 
and Rfiat = 0.03 pc. The dashed blue curve shows a Gaussian 
fit to the central part of the profile with a FWHM = 0.07 pc, 
which is considered as the inner width of the filament. The dot- 
dashed horizontal line indicates the value of the average column 
density along the spine while the dotted line indicates the average 
filament column density. 


clouds, namely the deconvolved width, the projected length 
(defined as the length of the longer spine build with the MST 
in the filament) not corrected for the inclination angle and 
the average column density. The median value of the dis¬ 
tributions are reported in Table The median value of the 
filament width is slightly lower, but still consistent with, the 
typical value of 0.1 pc found in a wide sample of filaments 


in other Gould Belt molecular clouds ( Arzoumanian et al. 
|2011[ ). Filaments are generally short with a median length 
of ~ 0.4 pc and only a few structures longer than 1 pc. The 
aspect ratio, however, is quite high, within the range from 
3 to 40 with median value of 7, 5 and 6 for Lupus 1, 3 and 
4, respectively, indicating the definite filamentary nature of 
the identified structures. The average filament column den¬ 
sity ranges from 6x10^° cm“^ to 7x10^^ cm~^. Noticeably, 
these values are in the lower part of the distribution of fila¬ 


ments column densities found in other molecular clouds (Ar- 
|zoumanian et aL]|2013p , and they are closer to the column 
densities found in ICS 146 and in the non-star-forming Po¬ 
laris region rather than to values typical of more active and 
massive regions such as Aquila in which more massive stars 
are formed. The derived average filaments column densities 
are affected by several sources of uncertainty: i) the < 20 
per cent error of the absolute flux calibration of the original 
photometric maps; ii) the uncertainty related to the assump¬ 
tion in the opacity law used to produce the column density 
maps that ranges from a factor of 1.5 for line of sights to¬ 
ward dense compact cores up to a factor of 2 for low column 
density regions (see Sect. 2.31; in) the unknown inclination 
angle of the filaments that however should be low for fila¬ 
ments having a high aspect ratio, as is the case the ones we 
have observed in Lupus; and iv) possible contamination by 
the background of a more diffuse medium in which filaments 
could be embedded and that has not been subtracted. In the 
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Lupus clouds, however, the diffuse background is negligible 
as can be seen in Fig.[^ indeed practically all the dense ma¬ 
terial at visual extinctions larger than 2 mag is arranged in 
filaments with part of filaments lying in regions of very low 
column density, at < 2 mag, especially in Lupus 1 and 
4 (see Figs. and 1^. Considering the low contribution of 

the last two points, we estimate an uncertainty on the de¬ 
rived average filaments column density of at most a factor 
of 2 . 

In summary, comparing the filaments in the Lupus com¬ 


plex with that of other molecular clouds | Arzoumanian et al. 
2013), and of the Galactic plane ( Schisano et al.||2014 ), the 
Lupus filaments have the typical width but a shorter length 
and, most importantly, a low column density. Approximat¬ 
ing the shape of filaments with a cylinder, a fixed observed 
width and a small column density imply also a small average 
volume density of the material arranged in filaments. In Lu¬ 
pus the low column density regime is not only valid for the 
material structured in filaments but is a property also of the 
diffuse material as testified by the probability distribution 
function analysis (see Sect. for a further discussion). 

It is worth noting that, in the three clouds both the 
range and the median values of the analysed filament param¬ 
eters are very similar (see Fig.[^. This result is in agreement 
with the observed similarity of the overall star formation 
properties of the three subregions of the Lupus complex. 
The general uniformity allows us to sum together the fila¬ 
ment samples from the three regions, in order to improve 
the statistical relevance of the subsequent analysis. 

3.2 Filament kinematics 

Information on the kinematics of filaments can be derived 
from spectroscopic observations of molecular lines. A good 
tracer of dense molecular gas is the CS molecule (critical 
density ~10® cm“®). Accordingly, we used the CS (2-1) line 
maps of Lupus 1 and 3 for a kinematical analysis of fila¬ 
ments; we do not have data for Lupus 4. The CS (2-1) line 
maps cover only the denser part of the Lupus 1 and 3 clouds, 
approximately the area above the level of 4 mag of visual ex¬ 
tinction. In the overlapping region the line emission, tracing 
the molecular gas, shows a morphological distribution which 
is very similar to the distribution of the dust inferred from 
the Herschel column density map. In fact, the filamentary 
structures identified in the Herschel maps are also present 
in the CS (2-1) line maps, see Fig. as example. 

The CS (2-1) channel maps of Lupus 1 and 3 (see fig. 
3 in Benedettini et al. (2012) and Fig. show that a sin¬ 
gle filament or a branch of a filament in the column density 
maps corresponds to a coherent kinematical structure in the 
channel maps. However, filaments of one cloud are not all at 
the same systemic velocity but they are at slightly different 
uisr, for example at 4.2 km s“^ and 4.7 km s“^ in Lupus 3 


and at 4.6 km s“^, 5.1 km s“^, and 5.8 km s“^ in Lupus 
1 (see Figs. 10 and 111. These differences, well above the 
spectral resolutionof the data of 0.1 km s“^, can be due to 
different inclination angles of the filaments or/and to small 
differences in the intrinsic velocity of the material compos¬ 
ing the filaments, possibly derived from the turbulence of 
the original molecular cloud. The kinematical coherence of 
filaments and the spread of 1-2 km s“^ of the velocity cen¬ 
troid among different filaments of the same cloud have been 


„-i 




Average filament column density (cm^) 

Figure 8. Distribution of the deconvolved width of the filaments 
(a), the length of the spine of the filaments (b) and the average 
filament column densities (c) of the three Lupus regions. The 
spatial resolution limit of the images is 0.026 pc for Lupus 1 and 
4 and 0.035 pc for Lupus 3. 


observed also in other low-mass star forming regions (e.g. 
Hacar fc Tafalla|2011 Hacar et al.|2013 l. 

From the analysis of the line profile we found that to¬ 
ward the central bright region of Lupus 3, the CS (2-1) line 
presents two peaks. This can be due to a self-absorption 
or to the overlap of two different gas components. We do 
not have the C^'^S ( 2 - 1 ) spectrum to verify which is the 
true case, however the fact that the peak of lines of differ¬ 
ent species, namely HC 3 N (10-9) and N 2 H"'" (1-0) with very 


similar critical density (~ 6 xl 0 ® cm ®) is coincident with 


only one of the CS peak (Benedettini et al. ][ 20 T 2 l ) supports 
the second hypothesis. In Fig. mr we show the map of the 
line centre of the two components Gaussian fitting: the up¬ 
per panel shows the distribution of the component at higher 
velocity (~ 4.7 km s“^), the lower panel shows the distri¬ 
bution of the component at lower velocity (~ 4.2 km s“^), 
highlighting two coherent structures of gas that superimpose 
each other in the central part of the map. A similar struc¬ 
ture has also been observed in the map of C ^*0 ( 2 - 1 ) line 
obtained with APEX (Benedettini et al., in preparation). 

No evidence of velocity gradients in the longitudinal or 
radial directions of the filaments is revealed in our spectro¬ 
scopic data, indicating that any possible filament contrac¬ 
tion, if present, must produce a velocity gradient, projected 
along the line of sight, below the 0.1 km s“^ spectral resolu- 
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Filaments in the Lupus molecular clouds 9 


Table 2. Median properties of the filaments distributions in the three Lupus clouds. 


Median deconvolved width (pc) 
Median length of the main spine (pc) 
Median average column density (10^^ 
Median aspect ratio 


Lupus 1 

Lupus 3 

Lupus 4 

0.05 

0.09 

0.07 

0.37 

0.31 

0.45 

1.2 

1.9 

1.6 

7.4 

4.9 

5.9 



i6^o'^oo"o og'^oo'o 

R.A.(J2000) 


Figure 9. CS (2-1) channel map of Lupus 3; first contour and contour steps are 0.1 K km s 
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Figure 10. Map of the central velocity of the Gaussian fitting of 
the CS (2-1) line in Lupus 1. The black lines indicate the spine 
of the filaments identified in the column density map. 
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Figure 11. Maps of the line centre derived from the Gaussian 
fitting of the CS (2—1) transition with two components in Lupus 
3: the upper panel shows the distribution of the component at 
higher velocity (~ 4.7 km s”^), the lower panel shows the distri¬ 
bution of the component at lower velocity (~ 4.2 km s“^). In the 
central part of the map the two components coexist. The spectral 
resolution of the data is 0.1 km s“^. The offset 0,0 corresponds 
to R.A.(J2000) = 16*’09™34?76 and Dec(J2000) = -39°06'54'.'l. 


tion of our data. Given the low value of this upper limit, we 
can confidently say that there is no evidence of contraction 
in the Lupus filaments at the present epoch. 

3.3 Filament arrangements 

Looking at the arrangement of the filaments in the Lupus 
clouds (see Figs. and [^, it appears that they are pre¬ 
dominantly arranged in the directions of the very long (a few 
parsec scale) principal filamentary structures. To verify this 
visual impression, we estimated the direction of the branches 
of the filaments as the position angle of the linear interpo¬ 


lation of the pixels located on their central spines. The dis¬ 
tribution of the position angles of the branches, measured 
counterclockwise from the west direction (north of west), 
for the three clouds is shown in Fig.[^ The distributions in 
each Lupus clouds are bimodal with two orthogonal peaks, 
one around 50° and the other around 140°. The bimodal 
distribution is less evident in Lupus 3 which can be due to 
the poor statistics of this cloud. The similarity of the po¬ 
sition angle distributions among the three clouds suggests 
that a common, probably external, event acted in shaping 
the material of the clouds and the formation of filaments. 
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Branch direction (degree) 


Figure 12. Distributions of the position angles (measured North 
from West) of the branches of the filaments of the three Lupus 
regions. The vertical bars indicate the orientation of the ambient 
magnetic field in each cloud, measured form the polarisation of 
the light from field stars. 


This common event could be the supernova explosion in the 
Upper-Sco sub-group of the nearby Scorpius-Centaurus OB 
2 association ( Humphreys|1978 1 that has been suggested to 
be the event that triggered the recent star formation in the 
Lupus complex ( Tachihara et al.|2001 |. 

Interestingly, one of the peaks in the position angle dis¬ 
tribution is similar to the orientations of the magnetic field 
in the regions (estimated from polarimetric observations of 
the field stars) that are: 124°-140°, 92° and 115° for Lupus 
1, 3 and 4, respectively ([Myers &: Goodman||1991| [Rizzo et| 
al. 1998| Matthews et al. 2014 b Assuming that the mag¬ 


netic field inside the cloud maintains a similar direction of 
that of the ambient field, and considering the large disper¬ 
sion associated to the mean value of the orientation of the 
magnetic field, we can assert that filaments in Lupus are 
mainly arranged along the direction of the magnetic field 
or in the orthogonal direction. This similarity is an indica¬ 
tion that the magnetic field may have also had a role in 
the filaments formation and may be important to determine 
their support, as also found in other star-forming regions 
(e.g., Palmeirim et al.]|2013 Busquet et al.|[2013 l. Recently 
[Li et al.| ( |2013| ) investigated the relation between the mag¬ 
netic field and the ISM distribution, and found that the few 
parsec length filamentary structures present in some of the 
Gould Belt molecular clouds preferentially lay either paral¬ 
lel, or orthogonal to the field direction. Our result in Lupus 
suggests that this is not only valid for the long filamentary 
structures, from a few to 10 pc scale investigated by [Li etj 


al. (20131, but also at the sub-parsec scale probed by our 


filaments that have length smaller than 1 pc. 


3.4 Association with compact sources 


Filaments have been demonstrated to be the preferred place 
for the formation of new stars (e.g. Andre et al.|2010 Poly- 


chroni et al.|2013 Konyves et al.|2015 1. However the under¬ 

lying physical processes that describe how the diffuse gas 
in the filaments fragments, and condenses into dense cores. 


remain unclear. The far-infrared Herschel maps are a pow¬ 
erful dataset for this kind of study since in these maps both 
filaments and compact dense cores from which stars form 
can be easily identified. In particular, the Herschel compact 
sources represent dusty dense cores that can have different 
natures, indeed they can be young stars still embedded in 
their parental envelope, cores without an internal source but 
dense enough to eventually form a star, or unbound cores 
that may eventually dissipate. The presence of a protostar 
that internally accretes from the dense core can be revealed 
by emission at wavelengths shorter than 100 pm that devi¬ 
ates from the grey-body spectral energy distribution typical 
of the emission at longer wavelengths. Therefore, a detection 
in the PACS 70 pm map can be used to identify the proto¬ 
stars. The further distinction among the starless cores of the 
gravitationally bound structures that are potential sites of 
star formation requires the estimate of the core mass. How¬ 
ever, a detailed catalogue of Herschel compact sources is out 
of the scope of the present paper and it will be the matter 
of a forthcoming paper (Benedettini et al., in preparation). 
Here, we refer to the preliminary list of Herschel compact 
sources published by Rygl et al. (20131. 

Looking at the distribution of those sources (see Fig. 
A.l of Rygl et al. [2013 \ it is clear that the majority of 
the starless cores (both gravitationally bound and unbound 
cores) are associated to filaments, as observed in other star 
forming regions (e.g. [Andre et ah 2010 Polychroni et al. 


2013[ Konyves et al.pOlSf In particular, cores are usually 


close to the central higher column density regions of the fila¬ 
ment (within one beam from the spine). Remarkably, in Lu¬ 
pus we find also filaments, or branches of long filaments, that 
do not contain compact sources. Filaments without cores are 
also observed in Polaris, a quiescent non-star-forming re¬ 
gion, while they are uncommon in more active and massive 
regions as Orion L1641. In Lupus those filaments without 
dense cores are mostly located in the lower column density 
regions of the cloud, indicating that the density fluctuations 
produced by the filament fragmentation may survive and 
grow preferentially in higher column density filaments. 


4 PROBABILITY DISTRIBUTION 
FUNCTIONS 

The probability distribution functions of column density 
(PDFs) for Lupus 1, 3, and 4 are displayed in Fig. [13[ Their 
shape is best fitted by a log-normal distribution at low col¬ 
umn densities and a power law at higher column densities, 
confirming earlier studies of PDFs in molecular clouds from 
extinction maps and Herschel column density maps (e.g. 


Kainulainen et aL][2009[ [Schneider et al.||2013| [Tremblin et[ 


al.[[2014 1 . However, a single log-normal fit for the very low 

column density range shows a deficit in the observed PDF 
which may account for a more complicated PDF structure. 
In particular, noise, definition of boundaries, image clipping, 
and projection effects make this low column density end dif¬ 


ficult to constrain by observations (Schneider et al. 2015 


Lombardi et al. 20151. In Table 3 we report the parame¬ 


ters of the best fitting for the three Lupus clouds. Both the 
widths of the log-normal part and the slopes of the power- 
law tail are in agreement with what is found in other Her¬ 


schel studies of low- and high-mass star forming clouds (Hill 
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Filaments in the Lupus molecular clouds 11 


Table 3. Parameters of the PDF fitting. For the log-normal func¬ 
tion; the column density at the peak <N>, the width cr and the 
maximum M; for the power law tail: the exponent s. 


<lV(H2)>(cm-2) 

Lupus 1 
5.1x102° 

Lupus 3 
9.7 x102° 

Lupus 4 
8.2 x102° 

a 

0.23 

0.17 

0.14 

M 

0.16 

0.23 

0.26 

s 

-2.4 

-2.4 

-2.6 


et al.|2011[ ISchneider et al.|2013] 2015; [A lves de Oliveira et| 


al.||2014f Tremblin et al.|2014| Konyves et al.||2015h 

Comparing the PDF of the three clouds, it is evident 
that the PDF of Lupus 1 extends to much lower column 
density, and has a peak at a value lower of about a factor 
of two compared to the other two clouds. This characteris¬ 
tic is partially due to the fact that in Lupus 1 we mapped 
not only the denser part of the cloud but also a large area 
of its more diffuse surrounding medium. Indeed, about the 
75 per cent of the field is composed of pixels at very low 
column density (Av < 1 mag), while it is not the case for 
the other two clouds. The very low column density of the 
diffuse material around Lupus 1 can be related to its pecu¬ 
liar position. This cloud is at the edge of an HI shell left 
over from a supernova event in Upper-Scorpius. The dy¬ 
namics of HI and CO toward Lupus 1 are consistent with 
the molecular cloud being compressed by the HI shell (de 


Geus 1992 Tothill et al. 20091. This scenario is supported 


by the morphology of the dense gas shown by the Herschel 
column density map, with the gas at > 1 mag confined 
in a thin layer parallel to the edge of the HI shell (see Fig. 

surrounded by the very low density medium swept by the 
wind of the supernova. On the other hand. Lupus 3 and 4 
are further from the ridge of the shell and their dense mate¬ 
rial is more distributed. The cloud morphology revealed by 
the Herschel column density maps once again indicates the 
role of the Scorpius-Centaurus OB association in the star 
formation history of the Lupus complex. 

A peculiarity of the PDFs of the Lupus clouds is that 
the values of the peak and the maximum of the column den¬ 
sities are systematically lower than those found in active 
regions such as Orion and Aquila, and appear to be much 
more similar to the values found in Chamaeleon and Po¬ 


laris (Schneider et al.]|2013 Schneider et al.||2015 Alves de 


Oliveira et al.|2()14 |, regions with low or null star formation 

activity. Remarkably, in the Lupus clouds the peak of the 
PDFs is at H 2 column density equivalent to visual extinc¬ 
tion less than 1 mag, the lowest value found so far in the 
PDF analysis of column density maps of star forming regions 
based on Herschel data, that usually peak at visual extinc¬ 
tion larger than 1.5 mag. These low values of the H 2 column 
density is an indication that in Lupus a part of the Hydrogen 
could be still atomic rather than molecular. Indeed, in Lupus 
the ratio of the mean column density of atomic Hydrogen 
( de Geus||1992 | over the peak the molecular Hydrogen dis¬ 
tribution is larger than one (~ 1.6 and 1.3 in Lupus 1 and 3, 
respectively) while, for example, in Orion B it is lower ~ 0.6 


(ISchneider et al.|2013 Ghromey, Elmegreen and Elmegreen 

198 ^ 

We separated the PDFs into ones with pixels on fila¬ 


ments and one with pixels off filaments (Fig. 131. In this 



1019 io2» 10=' 10 ^ 10^' 




f .I I ■ 'I'-.I_^.I: 'I . 

I019 1029 10=' 10^^ 10^' 

N(H2) Column Density [cm 


Figure 13. PDF of the Lupus 1, 3 and 4 clouds. The black his¬ 
togram represents the distribution of all the pixels of the map, the 
red, dashed histogram the distribution of the pixels belonging to 
filaments excluding pixels around dense cores and the blue dotted 
histogram the distribution of the pixels not belonging to filaments 
excluding pixels around dense cores. The log-normal and power 
law fit to the PDFs are also shown as continuous black lines. 


analysis we excluded pixels specifically belonging to dense 
cores, masking all the pixels within the HPBW from their 
reference positions. Such exclusion allows us to separate the 
contribution of the cores from that of the filaments at the 
high column density end of the distribution, where we ex¬ 
pect a strong contribution from gravitationally bound struc¬ 
tures. The off-filament PDF shows a log-normal distribution 
at low column densities but still has a power-law tail at 
higher column densities, albeit steeper than the total PDF. 
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12 M. Benedettini et al. 


On the other hand, on-filament PDF makes up the bluk of 
the power-law tail at higher column densities, confirming 
that the higher density medium is almost confined in fila¬ 
ments. A similar result has also been found by |Schisano et| 
al. (20141 for filaments in a portion of the Galactic Plane 


and by Konyves et al. (20151 in Acquila. Since deviations 


from the log-normal distribution in the form of a power law 
at high column densities are predicted for self-gravitating 


clouds ( Klessen]2000 Ballesteros-Paredes et al.|2011| , what 
we observe is an indication that self-gravity is important in 
some non-cores portions of the filaments and that locally 
the denser material contained in filaments, but still not in 
confined in a core, may undergo a gravitational collapse. 
This hypothesis can be also supported by the morphology of 
the material that composes the power-law tail of the PDFs. 
These pixels can be spatially identified by plotting the con¬ 
tour of the column density at which the on-filament dis¬ 
tribution deviates from the off-filaments one. This column 
density threshold is around 7x10^^ cm~^. The structures 
showed by this contour, although sitted on filaments, are 


rather clumpy (see Fig. 141 indicating that they are likely 


affected by gravitational collapse. A further indication of 
gravitational collapse in the identified structures is the link 
between their radial profiles and the index of the power-law 
in the PDF. A filamentary collapsing structure would give 
a radial profile oc r~°‘ where a is related to the exponen¬ 
tial of the PDF power law tail s through the relation a=l- 
1/s. Assuming for s the fitted values given in Table the 
corresponding filament profile would be very shallow with 
0=1.38, 1.40 and 1.42 for Lupus 1, 3 and 4, respectively. 
By contrast, clumpy collapsing regions would have a radial 
profile r““ with o= 1-2/s, namely o=1.77, 1.80 and 1.83 for 
Lupus 1, 3 and 4, respectively. These profiles fit very well 
with the classical values given from gravitationally-bound 
(free-fall and collapsing) structures with o between 1.5 and 
2. Furthermore, we checked the radial profiles of the regions 
above the 7x10^^ cm“^ threshold in Lupus 1 and 3 (an ex¬ 
ample is shown in Fig. 141 and we found < o >=1.89 for 


Lupus 1 and < o >=1.95 for Lupus 3, indicating local and 
clumpy collapsing regions rather than filamentary regions. 


5 DISCUSSION 

5.1 Lupus: a prototype of low column density star 
forming region 

Analysing the morphological distribution of the dense ma¬ 
terial in the Lupus complex, we find that almost all the ma¬ 
terial at visual extinctions larger than 2 mag is arranged in 
filaments and that some filamentary structures are present 
also in less dense region, e.g. in zones at Av < 1 mag. The 
Lupus clouds are characterised by having low column densi¬ 
ties, both in the dense material arranged in filaments, whose 
average column density is ~1.5xl0^^ cm“^, as well as in the 
diffuse medium, as seen from the peaks in the PDF distri¬ 
butions between 5xl0^°cm“^ and 10^^cm“^. The absence 
of very high column density gas and the predominance of 
low column density gas of course have an impact on the 
typical mass of the stars formed in Lupus that is quite low. 
Lupus is in fact one of the main low-mass star-forming com¬ 
plexes and its stellar population is dominated by mid M-type 


Lupus 3 



0.0 0.2 0.4 0.6 0.8 I.O 1.2 1.4 1.6 
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Figure 14. Top: A zoom of the column density map of the cen¬ 
tral region of Lupus 3. The contour of 7x10^^ cm~^ is drawn, 
showing that the higher column density material is confined in 
rather clumpy structures. Bottom: radial profile of clump 1. We 
fitted the external part with a function of the form N(H 2 ) oc r““, 
obtaining a best-fit value of a= 1.9 indicative of a gravitationally 
bound core. 


stars (Hughes et al. 19941 with a relative fraction of low- 
mass stars greater than other low-mass star-forming regions 
such as Taurus. In comparison to other molecular clouds, the 
star formation regime in Lupus, in terms of star formation 
rate and stellar clustering, represents an intermediate case 
between heavily clustered, active and more isolated, quies¬ 
cent star forming regions. Our study reveals that also the 
diffuse matter from which the stars are formed shows the 
same property: both the PDFs and the distribution of the 
filaments column density in Lupus are more similar to the 
values found in the quiescent Polaris region, rather than in 
the more active and massive Aquila and Orion star-forming 
regions. 


5.2 Can thermally subcritical filaments contain 
presteliar cores? 


A key physical parameter for filaments is the mass per unit 
length. Theoretical models of isothermal infinite cylindrical 
filaments confined by the external pressure of the ambient 
medium, predict the existence of a maximum equilibrium 
value to the mass per unit length ( Ostriker|1964 1 : 

2kT 


Afline,max — 


fj.mnG 


= 16.4 


10 K 


M© pc 


(3) 


where n is the mean molecular weight, mu is the mass of the 
hydrogen atom, k the Boltzmann constant, and G the grav- 
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Branch mass per unit length (Msun pc"') 

Figure 15. Distribution of the average mass per unit length of 
the filament branches.The vertical bar indicates the critical value 
of 16 Mq pc“^. 



k (|xm) 

Figure 16. SED of the Ll-11 core and grey-body fitting with T 
= 10 K and /3=2. The error associated to the fluxes corresponds 
to 20 per cent. See Fig. A.2 in |Rygl et al.| ( |2013^ for the five 
Herschel maps of the Ll-11 core. 


itational constant. When the filament mass per nnit length 
is higher than this maximnm value the filament is not ther¬ 
mally supported and it will collapse. In Fig. |15| we show the 
distribution of the masses per unit length of the branches of 
the filaments in Lupus. We found that almost all the filament 
branches have mass per unit length lower than the critical 
value of 16.4 Mq pc“^ for radial gravitational collapse and 
fragmentation of a cylindrical, isothermal filament with a 


typical temperature of 10 K. Previous studies (Andre et al. 


2010 2014; Arzoumanian et al. 20131 have suggested that 


only filaments with mass per unit length above this critical 
value (thermally supercritical filaments) are associated to 
star formation activity. However, in Lupus there is an evi¬ 
dent star formation activity even if the large majority of fila¬ 
ments are subcritical. More specifically, we have found that 
some subcritical filaments contain starless cores that have 
been classified as gravitationally bound cores by |Rygl et al.| 
(20131, since they satisfy the condition Mcore ^ Mbe, where 


Mbe is the Bonnor - Ebert mass of a core with equivalent 
size and temperature of the observed core. The gravitational 
bound cores are considered good prestellar core candidates 


in which new stars may form. In Rygl et al. (20131 the mass 


of the sources have been derived from a grey-body fitting 
of the SED of the fluxes measured the Herschel maps with 
the CUTEX algorithm ( Molinari et al.|2011 1. Anyway, disen¬ 
tangling the emission of a compact source from the emission 
of the embedding filament in the Herschel maps is not a 
trivial task. The most critical factors are the fit of the com¬ 
plex background and of the source size that introduces large 
uncertainty in the core mass estimate and therefore in es¬ 
tablishing if the core is gravitationally bound or not. For 
this reason it is of fundamental importance, when estimat¬ 
ing the flux of the sources in a complex region like molecular 
clouds, to compare the results of at least two completely in¬ 
dependent methods to have an idea of the robustness of the 
measured fluxes. 

To this aim we selected one well defined compact source 


in Lupus 1 - the prestellar core Ll-11 in Rygl et al. (20131 at 
R.A.(J2000) = 15*'42“5a6 and Dec(J2000) = -34n3^5'.'4 
- and we extracted the fluxes at 250 /rm, 350 /rm and 500 
/rm in our new Herschel maps by using the hyper algorithm 


( Traficante et al!]|2015 I which computes aperture photome¬ 
try instead of synthetic photometry like CUTEX. The corre¬ 
sponding SED is reported in Fig. |16| together with a grey- 
body fitting with fixed T = 10 K and (? = 2. The best fit 
core mass is 0.25 Mq, similar to the Bonnor - Ebert mass of 
a core with the same size of 36 arcsec and temperature of 10 
K, Mbe = 0.26 Mq. In conclusion, the hyper photometry 


confirms the result previously found by Rygl et al. (20131 


by using the CUTEX photometry that this source is a good 
candidate prestellar core. 

Interestingly, this candidate prestellar core is associated 
to a subcritical filament with a mass per unit length of 5 Mq 
pc”'^. Note that the filament can be still classified as sub¬ 
critical even if we consider the largest uncertainty quoted 
for the column density, namely a factor of 2. We stress that 
our estimate of the filament mass per unit length is a value 
averaged along the branch. In other words, we are assuming 
that the total mass measured in the entire structure was ini¬ 
tially uniformly distributed along the filament branch when 
it formed. Therefore, even if the branch is on average sub¬ 
critical, the presence of a gravitationally bound core is a 
indication that locally, in a smaller portion of the filament, 
the density is enough to become gravitationally unstable. 
The demonstration of if and how much the association of 
prestellar cores with thermally subcritical filaments is statis¬ 
tically significant requires the complete catalogue of prestel¬ 
lar cores in this region and it is deferred to a forthcoming 
paper (Benedettini et al., in preparation). 


5.3 The role of the magnetic field 


In Sect. 13.31 we showed that the filaments orientation in Lu¬ 
pus is correlated with the local magnetic field. The action 
of the magnetic field on filaments has been considered by 


models of Eiege & Pudritz (20001 and Eiege et al. (20041. 


In these models, the ratio between the mass per unit length 
and the virial mass per unit length Mune.vir = ‘iFlatlG is 
a key parameter of a filament that indicates whether it is 
in pressure equilibrium (Mune/Mune.vir ~1) or a poloidal 
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(Mline/Mline.vir >1) Or toroidal (Mline/Mline.vir <1) field 

is required to confine the filament. We calculated the virial 
mass per unit length by using the total velocity dispersion 
of the particles 


— \/crt\erm(f‘) + '^Lnth 
where the thermal dispersion is 


n'therm(/^) — 


kT 

firriH 


and the non-thermal velocity dispersion is given by 


<^nonth 






(4) 


(5) 


( 6 ) 


where pobs is the mean molecular weight of the observed 
molecule. 

We derived the non-thermal contribution for the fila¬ 
ments in Lupus from the line width of the CS (2-1) line 
(/roba=44), assuming that the line is optically thin. The line 
was observed with the Mopra telescope ( [Benedettini et SI] 
2012), at a spatial resolution of 35 arcsec, similar to that 


of our column density maps. The typical FWHM is about 
0.8 km s“^ and 0.4 km s“^ for Lupus 1 and 3 respectively 
that correspond to a total velocity dispersion of 0.39 and 
0.25 km s“^, values in agreement with the typical values 
found in thermically subcritical filaments ( [Arzoumanian et] 
al. [2013 ). The corresponding virial masses per unit length 


(70 Mq pc“^ for Lupus 1 and 30 Mq pc“^ for Lupus 3) 
are higher than average masses per unit length measured 
in Lupus indicating that the observed filaments are gravita¬ 
tionally unbound structures. Therefore, it is possible that an 
external agent will be required to provide confinement of the 
filaments. Since we find that Mune/Mune,viT <C 1 a binding, 
toroidal dominated, magnetic field is possibly confining fil¬ 
aments rather than just the external pressure, in agreement 


with what found also in other filamentary clouds (Fiege & 
Pudritz[[2000 ). 


How filaments form from the diffuse material of the 
molecular clouds is still unknown. Models for filament for¬ 


mation invoke shock compression (e.g., Heitsch & Hartmann 


2008 Padoan et al. 2001), shear and energy dissipation in 


large-scale MHD turbulence flow (e.g., [Hennebelle|2013P re¬ 
quiring the presence of strong magnetic fields. The first sce¬ 
nario foresees randomly distributed filaments without any 
correlation between the directions of clouds and mean mag¬ 
netic field directions (Padoan et al. 2001). In the second sce¬ 
nario, however, the competition between gravitational and 
turbulent pressures in a medium dominated by magnetic 
fields will shape the cloud to be elongated either parallel 
or perpendicular to the magnetic field when reaching equi¬ 
librium. In contrast to the shocked cloud model, stream¬ 
ing motions are not expected in the equilibrium stage of 
the magnetic field dominated model. In Lupus we observe 
a significant alignment of the filament directions with the 
ambient magnetic field and the absence of gas motions in¬ 
side the filament, both indications of a significant role of the 
local magnetic field in the formation of filaments. On the 
other hand, self-gravity seems not to play a major role in 
the filaments formation and fragmentation. Indeed, most of 
the filament branches in Lupus have average masses per unit 
length smaller than the critical value for radial gravitational 


collapse and no evidence of contraction has been seen in the 
gas kinematics. 


5.4 Kinematical structure 

The kinematical analysis of the filaments identified in col¬ 
umn density maps of the Lupus molecular clouds showed 
that they are composed by a network of short, coherent 
kinematical structures at slightly different velocity. Simi¬ 
larly, the (1-0) map of the B211/3 filament in Taurus 

revealed that the mostly continuous 10 pc long structure 
( Palmeirim et al.|2013| , visible in the column density map, 
is actually composed of many distinct components of gas, 
each of which has a coherent thermal velocity profile and 
typical length of 0.5 pc ( Hacar et al.[[2013 |. The inspection 
of the SPIRE 250 p.m map, where the small scale structure 
of the dust emission is better defined thanks to the higher 
spatial resolution and sensitivity, shows and excellent cor¬ 
respondence between the dust continuum B211/3 filament 


and the velocity-coherent fibres traced in C^“0 (Andre et 


al.|2014 l. The median length of the coherent kinematical fil 
aments observed in Lupus is similar to the typical length of 
the fibres observed in Taurus. 

This structure of interwined fibres was not observed in 
intermediate- and high-mass star-forming regions that have 
filaments usually thermally supercritical and with much 


higher column densities than observed in Lupus (Arzou- 
[manian et al.[[2013[ [Hill et 5(][201l[ ). Moreover, the kine 
matical analysis of some massive filaments have shown the 
presence of velocity gradients indicative of contraction and 
gravitational collapse ( [Schneider et al.[2010||[Kirk et al.[2013[ 
Peretto et al.[[2014 ), not observed in the Lupus low density 
filaments. 


6 CONCLUSIONS 

We have presented column density maps derived from the 
SED fitting of the fluxes at four far-infrared bands ob¬ 
served by Herschel for the Lupus 1, 3, and 4 star-forming 
regions. From these maps we extracted the filaments and 
calculated their physical properties. The Lupus clouds are 
characterised by low column densities both for the denser 
material assembled in filaments, whose average column den¬ 
sity is ~1.5xl0^^ cm“^, and from the diffuse medium, as 
testified from the peak of the PDFs between 5x10^° cm“^ 
and 10^^cm“^. This low column density regime of course re¬ 
flects upon on the mass of the formed stars in Lupus that 
indeed is low. 

Comparing our filaments catalogue with the prelimi¬ 
nary list of prestellar and protostellar sources identified by 


Rygl et al. (2013) in the Herschel maps, we find that most 


of the sources are located in filaments, in particular close to 
the central densest parts, testifying that filaments are frag¬ 
mented and contain the seeds for the formation of new stars. 
However, we find that, if we consider some average physical 
properties such as the average mass per unit length, most 
of the filaments in Lupus turn out to be gravitationally un¬ 
bound, having average mass per unit length lower than their 
respective virial critical masses. In addition, they seem to be 
stable self-gravitating structures, having average masses per 
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unit length lower than the maximum critical value for ra¬ 
dial gravitational collapse of an isothermal infinitely long 
hlament. Indeed, no evidence of Hlament contraction has 
been seen in the gas kinematics and each Hlament has an 
independent, coherent velocity structure. 

On the other hand, prestellar cores and protostars are 
observed in some of these supposedly globally subcritical 
filaments. This occurrence is an indication that in the low 
density regime the critical condition for the formation of 
stars can be reached only locally and it is not a global prop¬ 
erty of the filament. In fact, we have found evidences of the 
predominance of self-gravity in the highest column density 
regions of filaments in the fact that the power law tail of the 
observed PDFs are composed mostly by the pixels of the 
filaments and that these pixels belong to roundish clumpy 
structures with a radial profile compatible with clumpy col¬ 
lapsing regions. 

Multiple observational evidences indicate the key role 
that the magnetic field plays in the filament formation and 
confinement in this region: i) the masses per unit length 
are lower than the virial critical value, ii) one of the two 
peaks of the distribution of the position angle of the branches 
corresponds to the direction of the magnetic field, and iii) 
the absence of streaming motions observed in the filament 
kinematics. 
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